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Pairing of fermions is ubiquitous in nature and 
it is responsible for a large variety of fascinating 
phenomena like superconductivity, superfluidity 
of 3 He, the anomalous rotation of neutron stars, 
and the BEC-BCS crossover in strongly interact- 
ing Fermi gases. When confined to two dimen- 
sions, interacting many-body systems bear even 
more subtle effects 1 , many of which lack under- 
standing at a fundamental level. Most striking is 
the, yet unexplained, effect of high-temperature 
superconductivity in cuprates, which is inti- 
mately related to the two-dimensional geome- 
try of the crystal structure. In particular, the 
questions how many-body pairing is established 
at high temperature and whether it precedes 
superconductivity are crucial to be answered. 
Here, we report on the observation of pair- 
ing in a harmonically trapped two-dimensional 
atomic Fermi gas in the regime of strong coupling. 
We perform momentum-resolved photoemission 
spectroscopy 2,3 , analogous to ARPES in the solid 
state 4 , to measure the spectral function of the gas 
and we detect a many-body pairing gap above the 
superfluid transition temperature. Our observa- 
tions mark a significant step in the emulation of 
layered two-dimensional strongly correlated su- 
perconductors using ultracold atomic gases. 

One essential mechanism to establish collective quan- 
tum phases in Fermi systems, such as superconductivity 
or superfluidity, is fermionic pairing. In the case of weak 
attractive interactions between spin-up and spin-down 
fermions in three dimensions, Cooper pairs form on the 
surface of the Fermi sea, even though in vacuum no two- 
particle bound state exists. Pairing and the condensation 
of pairs go hand in hand. However, in this weak coupling 
limit the gap and hence the critical temperature T c are 
exponentially small. The fundamental and technological 
quest for increasing the superconducting transition tem- 
perature has led to the discovery of correlated materials 
which depart from weak coupling theory. These are com- 
plex compounds which combine strong interactions and 
large quantum fluctuations due to a two-dimensional ge- 
ometry, both of which yet hinder a precise understand- 
ing and thorough theoretical modeling, even at zero tem- 
perature. Strongly correlated, two-dimensional, layered 
superconductors have displayed a wealth of new and un- 
expected phenomena. Among those is the peculiar ob- 
servation of a suppression of low-energy weight in the 
single particle spectral function above the critical tem- 
perature and hence without direct link to superconduct- 
ing behavior—. The origin of this pseudogap phenomenon 



so far lacks full understanding and the implications of the 
pseudogap for the appearance of superconductivity have 
remained an open question. 

Ultracold atomic gases provide an ideal testing ground 
to uncover the physics of strongly correlated Fermi sys- 
tems because of the exceptional capability to tune the 
system parameters. Experiments in three-dimensional 
Fermi gases have observed fermionic many-body pair- 
ing in the crossover regime from BCS-like pairing to 
Bose-Einstein condensation (BEC) of dimers (for a re- 
view see£). Neither in the BCS- nor in the BEC-limit a 
pseudogap is expected in three dimensions. In between, 
i.e., for kp\a\ ^> 1, the gas is in the unitary regime of 
strongest interactions. The Fermi wave vector is kp and 
a is the s-wave scattering length between spin-up and 
spin-down particles. Whether a pairing pseudogap phase 
exists here or whether a Fermi-liquid state dominates is 
still debated^. 

When the motion of particles is constrained to a two- 
dimensional plane, a bound state with binding energy 
Eb exists for any value of the three-dimensional scat- 
tering lengthkikil. This is in contrast to the situation 
in three dimensions where weakly bound molecules start 
to appear in the unitary regime of strongest interaction. 
The presence of the two-body bound state in vacuum is a 
necessary and sufficient condition for many-body pairing 
in a medium in two dimensions— and a small binding 
energy Eb -C Ep corresponds to the attractive BCS-like 
regime with interaction parameter \n.(kpa,2D) > 1- Here, 
Ep denotes the Fermi energy, a^o = h/ \frnE~B > is the 
two-dimensional scattering length, m is the mass of the 
atoms and h is Planck's constant divided by 2ir. At zero 
temperature, BCS-type superfluidity has been predicted 
in the attractive two-dimensional mean-field regime^—. 

Qualitatively, we understand the behaviour at strong 
coupling and finite temperature by considering a complex 
order parameter- A(x) = |A(x)|e j6l ( x ). The phase 6*(x) 
may fluctuate as a function of the spatial coordinate x, 
for example driven by thermal fluctuations and quantum 
fluctuations enhanced in low-dimensional systems. For 
temperatures below the Berezinskii-Kostcrlitz-Thouless 
transition temperature Tbkt — T c , the superfluid gap 
A sc = (A) is larger than zero. For higher temperature 
T > Tbkt phase fluctuations destroy long range order 
(i.e. (e ie(x) ) = 0), however, the modulus A PG = (|A|) 
remains finite up to the pairing crossover temperature 
T* , which is on the order of Eb- The regime between 
Tbkt and T* is referred to as the pseudogap regime and 
is still little understood. The transition between the nor- 
mal phase and the pseudogap phase is not associated 
with breaking a continuous symmetry but is a crossover 
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FIG. 1. Measuring the spectral function, a: Measured pho- 
toemission signal at ln(fcF<i2D) = 0. The image is averaged 
over six repetitions and we plot log 10 of the optical density, 
b: Energy distribution curves for various values of the mo- 
mentum hk. The data of the pairing peak (lower feature) 
are fitted with a modified Gumbel distribution and the peak 
positions (red crosses) indicate the backbending of the dis- 
persion curve. The data of the atomic peak are fitted with a 
Gaussian. 



phenomenon. This scenario is different from standard 
weak-coupling mean-field BCS theory in which the su- 
perfluid gap is destroyed by pair breaking and T c and 
T* are approximately equal. The temperature range for 
accessing the pseudogap phase, i.e., the difference be- 
tween T* and T c , is particularly large in two dimensions 
since enhanced quantum fluctuations and second-order 
interaction effects^ suppress Tbkt with respect to the 
mean- field result. The pseudogap regime is connected by 
a crossover to a Bose-liquid regime of local pairs when 
ln(fci?a2r>) crosses zero or, equivalently, Eb — IE p. 

Here, we study two-dimensional Fermi gases for T > 
Tbkt in the strongly interacting regime | ln(fci?a2£))| < 1, 
in which the size of the pairs is comparable to the mean 
interparticle spacing. Hence, Pauli blocking and other 
many-body effects act on the pairing and a key question 
is, whether a pseudogap phase exists above the super- 
fluid transition temperature^^! T B kt ~ O.lTp. We in- 
vestigate many-body pairing using momentum-resolved 
photoemission spectroscopy, see methods. To this end, 
a long- wavelength photon of energy Ml creates a single 
particle excitation in the two-dimensional gas, which we 
detect at the momentum k. This measures the spec- 
tral function A(Q, k) multiplied by the fermionic oc- 
cupation function. This technique was originally de- 
veloped for the study of solid state materials, such as 
cuprate superconductors^, but recently has been suc- 
cessfully adapted to cold atomic systems^. Figure la 
shows a measured spectral function as a typical dataset 
used in this paper. The spectrum exhibits two features: 
The upper branch near zero energy corresponds to un- 
paired atoms, most likely from low-density regions of the 
trap with vanishingly small interaction between the free 



atoms2£. The lower branch corresponds to the pairing 
signal with a non-trivial dispersion in the spectral func- 
tion. In figure lb we plot the spectra A{£1) for constant 
values of fc, which show clearly asymmetric line shapes 
as expected for pair breaking into a continuum, together 
with our fits (see Methods). 

We analyze the spectral function for various interac- 
tion strengths and temperatures. In the attractively 
interacting regime ln(kpa2D) > the size a,2D of the 
pairs is larger than the mean interparticle spacing kZ . 
Therefore, at low temperature no isolated dimers exist 
in the sample but a many-body pairing gap opens in the 
spectrum, conceptually similar to a BCS-paired state in 
three dimensions. The zero-temperature limit of the two- 
dimensional BEC-BCS crossover has been investigated in 
the framework of mean-field theorjii^. This has led to 
the remarkable prediction that at zero temperature the 
condensation energy per particle E = A 2 SC /2E F at k = 
should be equal to the two-body binding energy in vac- 
uum Eb for all interaction strengths. 

The expected rf spectrum can be qualitatively esti- 
mated based on mean-field theory of a BCS state at 
zero temperature. Many-body pairing causes an energy 
threshold E t h and this energy has to be supplied by the rf 
field to create a single-particle excitation. Ignoring final 
state interactions and assuming two-dimensional geome- 
try with a constant density of states, the momentum in- 
tegrated lineshape has the form r(f2) oc Q(Ml — E t h)/^ 2 ■ 
Here, 0(x) is the Heaviside function. We investigate 
many-body pairing at low temperatures by studying the 
energy threshold E t h of the pairing peak in the photoe- 
mission spectrum at momentum k — 0. This quantity 
is closely linked to Eq and its disappearance, for exam- 
ple as temperature is increased, indicates the transition 
into a normal phase. Investigating E t h at k = 0, rather 
than for example at kp, has the advantage that averaging 
over the inhomogeneous spread of kp and /i of the ap- 
proximately 30 two-dimensional gases has a small effect 
because the mean-field prediction at zero temperature is 
not explicitly dependent on fi and kp. At the tempera- 
ture T/Tp — 0.27 (see Figure 2a, top row), we observe 
a negligible contribution of free atoms and a sharp on- 
set of the pairing peak. As we increase the temperature 
of our sample (see Figure 2a, bottom row) the spectral 
weight at low energies increases because the pairing peak 
broadens and the threshold shifts. For \n(kpa,2D) = 0.8 
at T/Tp = 0.65, for example, the broadening creates a 
significant amount of low-energy spectral weight and the 
threshold of the spectrum even moves to positive ener- 
gies. The latter result shows that the energy gap has 
completely disappeared. 

In Figure 2b we show the measured temperature de- 
pendence of E t h for various interaction strengths on the 
BCS side of the resonance ln(kpa 2 D) > 0, where no lo- 
cal dimers exist in the sample. We observe that E t h 
moves towards zero with increasing temperature, which 
indicates that the many-body pairing gap vanishes. We 
compare our data to finite temperature mean-field theory 
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FIG. 2. Pairing in the strongly interacting two-dimensional Fermi, a: Energy distribution curves A(k — 0, fi) for different 
values of ln(fcFa2i)) and temperatures. For low temperatures hardly any unpaired atoms exist. The solid line is the fit described 
in the text, b: Energy threshold of the pairing peak as a function of temperature for different interaction strengths. The solid 
lines are predictions of mean-field theory^ for the quantity \A(T)\ 2 /2Ef in two dimensions without free parameters. The 
shaded area marks the uncertainty of our determination of \n(kFa,2D), which is ±0.1. The error bars show the l-cr fit errors. 
They become larger towards larger values of ln(fc_Fa2o) because of the division by the exponentially smaller binding energy Eb- 



in two dimensions^. Overall, the numerically calculated 
energy Eq — | A(T)\ 2 /2Ep reproduces our experimental 
data for E t h without free parameters (solid line in Figure 
2b). However, some discrepancies remain. For example, 
for the coldest temperatures the average of E t h / Eb over 
all interaction strengths is E t h/EB = (0.84 ± 0.05), ap- 
proximately 15% smaller than predicted by mean-field 
theory. The deviation could stem from beyond mean- 
field effects provoked by our two-dimensional geometry 
and interaction energy shifts. It could be more pre- 
cisely computed using Quantum-Monte Carlo calcula- 
tions, which, however, up to now exist only for zero 
temperature 2 ^. From our data we conclude that we have 
realized a many-body pairing phase in two dimensions 
and observed how it vanishes as the temperature is in- 
creased to approach the pairing crossover temperature 
T* . Our data are taken above the critical temperature for 
superfluidity Tbkt in two dimensions, which is predicted 
to be at T/Tp ~ 0.1 in the strong coupling regime. We 
have verified this experimentally by employing a rapid 
magnetic field sweep across the Feshbach resonance to 
project the possibly condensed pairs onto deeply bound 
molecules and then image the molecular condensate af- 
ter timc-of- flight 26 . While this procedure allows us to 
observe three-dimensional Fermi condensates, in our two- 
dimensional samples no condensation was detected. We 
therefore identify our phase as a pairing pseudogap phase. 

In addition to the energy threshold at zero momen- 
tum, we also investigate how the quasiparticle disper- 
sion is affected when changing temperature and interac- 
tion strength. We extract the dispersion data from the 
peak of the spectral function E pea k(k) = max$i [A{k, £1)] 
(see Figure lb). In the crossover region and on the 
BEC side of the resonance, the dispersion shows a pro- 
nounced back-bending, which we display in Figure 3a 
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FIG. 3. Interaction dependence of the quasiparticle disper- 
sion, a: Dispersion for different interaction strengths at 
T/Tf = 0.27. In the strongly interacting crossover regime, 
we observe back-bending of the dispersion relation which van- 
ishes in the regime of small interaction. The solid line is the 
fit of a BCS-like dispersion relation defined in the text, b: 
Wave vector k* at which the back-bending of the spectral 
function occurs, extracted from a fit to the dispersion rela- 
tion. The dashed line indicates the Fermi wave vector, c: 
Effective mass, extracted from a fit to the dispersion relation. 
The solid lines in b and c are the prediction of a thermal 
singlet model. The error bars reflect the fitting errors. 



for various interaction strengths. This shape is reminis- 
cent of the hole dispersion of a BCS state and qualita- 
tively similar to that of a three-dimensional Fermi gas 
near unitarity'. We fit the dispersion relation E(k) = 

H - J [h 2 /(2m*)(k 2 - k* 2 )] 2 + A 2 for k < 1.5k F using 
the zero temperature mean-field prediction for the chem- 
ical potential^ to extract the parameters effective mass 
to* and the wave vector k* , at which the backbending oc- 
curs. In Figures 3b and 3c we display the results of k* and 
to*, respectively, which show that k* remains similar to 
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kp = 8.1 /im _1 across the whole crossover. We note that 
the temperature of the cloud increases on the far BEC 
side, possibly due to heating introduced by the nearby 
confinement-induced p-wave resonance^ 7 -. We model our 
data on the BEC side of the resonance with a thermal 
ensemble of singlet pairs using the exact expression of 
the pair wave function in two dimensions and the ex- 
perimentally determined temperatures^ (see Methods). 
The approximation of phase-disordered pairs in the nor- 
mal state includes the correct short-range physics but ne- 
glects interactions between pairs as well as Pauli-blocking 
in the dissociation process. Our numerical data show 
that for experimentally relevant parameters the thermal 
singlet model can give rise to a back-bending feature of 
the photoemission spectrum with similar values of m* 
and k* as in the experiment (see solid lines in Figure 3b 
and c). A contribution to the back-bending effect from 
incoherent pairs has been predicted previously in three 
dimensions^. We note that, despite working reasonably 
well for k* and m*, the thermal singlet model makes in- 
correct predictions for the energy threshold in Figure 2. 
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FIG. 4. Temperature dependence of the quasiparticle disper- 
sion. a-c:Temperature dependence of the backbending sig- 
nature of the quasiparticle dispersion at \n(kF<i2D) = —0.2. 
d: Evolution of the pairing temperature Tg as a function of 
the interaction parameter \n(kFa.2D) ■ The dashed line is the 
theoretical prediction of T* from mean-field theory, scaled by 
a factor of 0.36. 



It is instructive to observe what happens to the disper- 
sion signal at elevated temperatures. To this end, we have 
performed experiments in the range 0.27 < T/Tp < 0.65 
in the strongly interacting crossover regime. We observe 
that the BCS-like dispersion relation converts into a free- 
particle like dispersion as the temperature is increased 
(see Figure 4a-c) and we determine the temperature T B 
at which this change takes place. Figure 4d shows the 
variation of T B (in units of Eb) in the crossover regime 
and on the BEC side. We find that the temperature T B 
is comparable to the pairing temperature T* of mean- 
field theory (dashed line in Figure 4d), however, a factor 
of 0.36 smaller. This discrepancy could suggest that the 
appearance of the back-bending feature in the spectral 
function, which has been interpreted as a signature for 
many-body pairing^, is mainly a qualitative evidence. 
In contrast, our spectroscopic measurements of the en- 
ergy threshold (see Figure 2) give quantitative results for 
many-body pairing and agreement with theory even in a 
regime where no back-bending feature is observable. 

In conclusion, we have observed a many-body pair- 
ing gap above the superfluid transition temperature of a 
strongly interacting two-dimensional Fermi gas. Our re- 
sults represent the first major step in emulating and un- 
derstanding pairing in two-dimensional, strongly corre- 
lated materials using cold atoms. In future experiments, 
an even closer cross-link to cuprate high-temperature su- 
perconductors could be provided by the inclusion of an 
optical lattice potential to realize the two-dimensional 
Fermi-Hubbard model. In this model, a pseudogap phase 
is expected for a slightly doped antiferromagnetic state 
and it is believed to precede exotic quantum states like 
the d-wave superconducting phase, which will be de- 
tectable by our technique. 



METHODS SUMMARY 

In our experimental setup 1 - 7 -, we prepare a quantum 
degenerate Fermi gas of 40 K atoms in a 50/ 50 mixture of 
the two lowest hyperfine states \F = 9/2, mp = —9/2) 
and \F = 9/2, mp = —7/2). We confine the quantum 
gas to two dimensions in a deep optical lattice formed 
by a standing wave laser field, preparing approximately 
30 layers. The interaction strength between spin-up and 
spin-down particles is tuned at a Feshbach resonance near 
202.1 Gauss. The photoemission measurement couples 
the \F = 9/2, mp = —7/2) state to the weakly inter- 
acting state \F — 9/2, mp — —5/2) using a radiofrc- 
quency photon of frequency SI with negligible momen- 
tum transfer. We measure the momentum distribution 
of the transferred atoms in a time-of-fiight experiment 
and average the absorption signal azimuthally to obtain 

A(k, ft), where k = Jk 2 x + kl. 



METHODS 

A. Experimental setup 

We evaporatively cool a 50/50 spin mixture of 40 K 
atoms in the \F = 9/2, m F = -9/2) = | - 9/2) and 
\F = 9/2, m F = -7/2) = | - 7/2) states of the hy- 
perfine ground state manifold 1 - 7 -. After reaching quan- 
tum degeneracy in a crossed-beam optical dipole trap 
with approximately 70000 atoms per spin state, we turn 
on an optical lattice potential in order to prepare two- 
dimensional Fermi gase a 17 ' 27 ! 29 ' 30 . The optical lattice 
is formed by a horizontally propagating, retro-reflected 
laser beam of wavelength A = 1064 nm, focussed to a 
waist of 140 fim. We increase the laser power over a 



5 



time of 200 ms to reach a final potential depth of up to 
Vi a t — 83 E rec , which is calibrated by intensity modula- 
tion spectroscopy. E rec = h 2 /(2mX 2 ) is the recoil energy. 
The trapping frequency along the strongly confined di- 
rection is uj = 2tt x 78.5 kHz. After loading the optical 
lattice, we adiabatically reduce the power of the optical 
dipole trap such that the atoms are confined only by the 
Gaussian intensity envelope of the lattice laser beams. 
The radial trapping frequency of the two-dimensional 
gases is ui± = 2tt x 127 Hz for Vi at = 83E rec and we 
confine on the order of 10 3 atoms per two-dimensional 
gas at the center of the trap. Along the axial direction 
we populate approximately 30 layers of the optical lattice 
potential with an inhomogeneous peak density distribu- 
tion. Approximately two thirds of the 2D layers with 
highest density dominate the measured signal and their 
relevant energy scales Ep, Eb, and A 2 /2Ep are more 
than an order of magnitude larger than the trapping fre- 
quency luj_. Therefore, finite particle number effects do 
not influence the measured signal. After evaporation, we 
adiabatically increase the interaction strength by lower- 
ing the magnetic field, at a rate of up to 0.25 G/ms, to 
a value near the Feshbach resonance at 202.1 G. We ap- 
ply a radio- frequency pulse near 47 MHz with a Gaussian 
amplitude envelope with a full width at half maximum 
of 230 /us to transfer atoms from the | — 7/2) state to the 
|F = 9/2, m F = -5/2) state. Atoms in the |9/2, —5/2) 
state have a two-body s-wave scattering length of 130 
Bohr radii with the | — 7/2) state and 250 Bohr radii 
with the | — 9/2) state 3 . We turn off the optical lattice 
100 fis after the radiofrequency pulse, switch off the mag- 
netic field and apply a magnetic field gradient to achieve 
spatial splitting of the three spin components in a Stcrn- 
Gerlach experiment. For each run, the magnetic field 
is calibrated using spin-rotation with an rf pulse of an 
imbalanced mixture on the | — 9/2)/| — 7/2) transition. 
The magnetic field accuracy deduced from these mea- 
surements is < 3mG. We measure the temperature by 
ballistic expansion of a weakly interacting gas, and the 
quoted numbers refer to the average of T/Tp across the 
whole sample. 



B. Determination of the energy threshold E t h of 
the energy distribution curve 

We fit our data with a double-peak fitting function 
comprising of a Gaussian for the atomic signal and a 
modified Gumbel function /(O) = aexp[— (f2 — ilo)/b — 
aexp(— (f2 — Qo)/(a&))] for the pairing peak. The param- 
eter fio measures the peak position and the parameters a 
and b measure skewness and width. For our further anal- 
ysis, we only use the peak position Qq, which does not de- 
pend on the line shape function used. From this fit we de- 
termine the maximum of the molecular peak v max = £Iq 
and the minimum between the atomic and the molecular 
peak v min . Between v x = v max and v 2 = v mm - 2 kHz 
we fit the data with a linear function and determine the 



zero-crossing of the linear extrapolation as the energy 
threshold E t h- We correct the obtained result for our 
spectral resolution of 1.5 kHz, obtained from the width of 
the Gaussian fits. The data are normalized to the two- 
body binding energy in vacuum which we obtain from 
the transcendental equation 6 

I la= I™ dU (l - ex P(~^W(M) \ 
Jo y/I^P \ v /(l-cxp(-2u))/(2u) J ' 

(1) 

Here, l z = y^h/mcu and a is the three-dimensional scat- 
tering length using the following parameters of the Fes- 
hbach resonance: Bo = 202.1 G, AB — 7G and a^G = 
174 as where a# is the Bohr radius. 



C. Thermal singlet model 

We model our data on the BEC side of the resonance 
with a thermal ensemble of singlet pairsS.. The expression 
for the wave function of the bound state in two dimen- 
sions is ipB(r) — \/2/a2DKo(r/a2D), in which K a (x) is 
the modified Bessel function, and for the scattering state 
is ipq( r ) = Ja{l r ) — l -^^-H t ^'\qr), in which Jq(x) is the 

Bessel function of the first kind and Hq (x) is the Han- 
kel function of the first kindi£. /(g) is the scattering 
amplitude between the state | — 7/2) and the final state 
| — 5/2). We compute the momentum resolved rf spec- 
trum for the dissociation from the bound state to the 
scattering state, averaging over a thermal distribution 
of the center-of-mass momenta of the initial pairs using 
Monte-Carlo sampling. From the momentum-resolved rf 
spectrum we calculate the effective mass m* and the wave 
vector k* using the same fitting routines as for the exper- 
imental data. This model of tightly bound pairs in the 
normal state includes the correct short-range physics but 
neglects many-body pairing, interactions between atoms 
and between pairs, as well as quantum statistical effects. 
Therefore, we do not expect quantitative agreement in 
the strongly interacting regime or on the BCS side of the 
resonance. 
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